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I discuss recent progress in understanding the high-energy evolution in QCD, which points 
towards a remarkable correspondence with the reaction-difltusion problem of statistical physics. 
This enables us to determine the asymptotic behaviour of the scattering amplitudes in QCD. 



1 Introduction 

Over the last year, an intense activity in the field of high-energy QCD has been triggered by 
the following observations: (i) the gluon number fluctuations in the dilute regime at low energy 
play an important role in the evolution towards gluon saturation and the unitarity limit with 
increasing energy^''^, (ii) the QCD evolution in the presence of fluctuations and saturation is a 
classical stochastic process which is similar to the 'reaction-diffusion' problem widely studied in 
the context of statistical physics and (iii) the relevant fluctuations are however missed^ by 
the existing approaches to non-linear evolution in QCD at high energy, namely, the Balitsky- 
JIMWLK equations These observations, together with their consequences, have entailed 
important conceptual clarifications and stimulated new ideas and theoretical constructions. 

The correspondence between high-energy QCD and statistical physics was in fact anticipated 
by the probabilistic structure inherent in previous approaches like the color dipole picture^^~^'^ 
and the color glass condensate ^' (the QCD effective theories at low and high gluon density, 
respectively). The recent developments in Refs.'^"^ made this correspondence more precise and 
also useful (in the sense of generating new results for QCD), first at the level of the mean 
field approximation — where the link ^ between the Balitsky-Kovchegov (BK) equation ^' ^'^ 
in QCD and the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation in statistical 
physics has shed a new light on the important phenomenon of geometric scaling ^"^ — , then in 
the analysis of the particle number fluctuations — where recent advances in statistical physics 
have enabled us to compute QCD scattering amplitudes under asymptotic conditions (very 
high energy and arbitrarily small coupling constant) These new approaches have elegantly 
confirmed and extended previous results obtained through direct studies in QCD ^' 

At the same time, it became clear that the correspondence with statistical physics cannot 
be used to also study the pre-asymptotic behaviour in QCD, that is, to compute scattering 
amplitudes for realistic values of the energy and the coupling constant. In that regime, which 
is the only one to be interesting for the phenomenology, one rather needs the actual evolution 
equations in QCD at high energy. As aforementioned, these equations should be more general 
— in the sense of also including the effects of gluon number fluctuations — than the previously 
known Balitsky-JIMWLK equations. So far, the relevant equations have been constructed^' 
only in the limit where the number of colors is large. An ambitious program, which aims at 



generalizing these equations to arbitrary values of Nc, is currently under way °. This effort 
led already to some important results — in particular, the recognition^"^' of a powerful 'self- 
duality' property of the high-energy evolution, and the construction of an effective Hamiltonian 
which is explicitly self-dual — , but the general problem is still under study, and the 
evolution equations for arbitrary Nc are not yet known. 

In my two succinct contributions to these Proceedings, I shall restrict myself to the large- A'c 
limit, which is quite intuitive in that it allows the use of a suggestive dipole language^' In this 
context, I shall rely on simple physical considerations to explain the correspondence between 
high-energy QCD and statistical physics (in this presentation), and then motivate the structure 
of the recently derived 'evolution equations with Pomeron loops' (in my other presentation^^). 

2 QCD evolution at high energy 

To put the theoretical developments into a specific physical context, let us consider 7* proton 
deep inelastic scattering (DIS) at high energy, or small Bjorken-x. We shall view this process in 
a special frame in which most of the total energy is carried by the proton, whose wavefunction 
is therefore highly evolved, while the virtual photon has just enough energy to dissociate long 
before splitting into a quark-antiquark pair in a colorless state (a 'color dipole'), which then 
scatters off the gluon distribution in the proton (see Fig. 1). The transverse size r of the dipole 
is controlled by the virtuality of 7* (roughly, ~ 1/Q'^), so for 3> ^qqd treat 
the dipole scattering in perturbation theory. But for sufficiently small x, even such a small 
dipole can see a high-density gluonic system, and thus undergo strong scattering. 

Specifically, the small-x gluons to which couple the projectile form a color glass conden- 
sate^' (CGC), i.e., a multigluonic state which is characterized by high quantum occupancy, 
of order l/ag, for transverse momenta below the saturation momentum Qs{x), but which 
becomes rapidly dilute when increasing k± above Qg- The saturation scale rises very fast with 
the energy^^, Q'^{x) x^^^ and is the fundamental scale in QCD at high energy. In particular, a 
small external dipole with size r <C l/Qs undergoes only weak scattering (since it couples to the 
dilute tail of the gluon distribution at large /c^), while a relatively large dipole with r ^ l/Qg 
'sees' the saturated gluons, and thus is strongly absorbed. 

In turn, the small-x gluons are produced through 'quantum evolution', i.e., through radiation 
from color sources (typically, other gluons) at larger values of x, whose internal dynamics is 
'frozen' by Lorentz time dilation. Let r = Inl/x denote the rapidity; it takes, roughly, a 
rapidity interval At ~ l/ccs to emit one small-x gluon; thus, in the high energy regime where 
asT 3> 1, the dipole meets with well developed gluon cascades, as shown in Fig. 1. Three types 
of processes can be distinguished in Fig. 1, which for more clarity are disentangled in Fig. 2. 

The first process. Fig. 2. a, represents one step in the standard BFKL evolution by 
iterating this step, one generates gluon ladders which are resummed in the solution to the 
BFKL equation^''. By itself, this mechanism entails a rapid growth of the gluon distribution 
(exponential in r), which however leads to conceptual difficulties at very high energy : (i) The 
BFKL estimate for the dipole scattering amplitude Tr(r) grows like a power of the energy, and 
thus eventually violates the unitarity bound r,-(r) < 1. (The upper limit Tt = 1 corresponds to 
the 'black disk' limit, in which the dipole is totally absorbed by the target.) (ii) The BFKL 
ladder is not protected from deviations towards the non-perturbative domain at low transverse 
momenta /c^ ^ ^qcd ('infrared diffusion'). With increasing energy, the BFKL solution receives 
larger and larger contributions from such soft intermediate gluons, and thus becomes unreliable. 

These 'small-x problems' of the BFKL evolution are both cured by gluon saturation^^ , the 
mechanism leading to the formation of a CGC: At sufficiently high energy, when the gluon 
density in the target becomes very large, the n — > 2 recombination processes illustrated in Fig. 
2.b start to be important and tame the growth of the gluon distribution. Such processes are 
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Figure 1: An instantaneous gluon conSguration in the proton wavefunction as 'seen' in DIS at small x. 
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Figure 2: Gluon processes which occur in one step of high energy evolution. 

included (to all orders) in the JIMWLK equation ''^^ , a non-linear and functional generalization 
of the BFKL equation which describes the evolution of the ensemble of gluon correlations in the 
approach towards saturation. (As manifest on Fig. 2.b, the standard 'gluon distribution', which 
is a 2-point function, gets coupled to the higher n-point functions via the recombination effects.) 
Remarkably, the JIMWLK equation is a (functional) Fokker-Planck equation, which describes 
the high-energy evolution as a classical stochastic process'^'^ — a random walk in the functional 
space of gluon fields. When applied to scattering amplitudes for simple external projectiles (like 
color dipoles, quadrupoles, etc.), the JIMWLK equation generates an infinite set of coupled 
evolution equations that were originally derived by Balitsky^. 

However, as recently noticed in Ref. ^, the Balitsky- JIMWLK hierarchy misses the 2 — > n 
splitting processes illustrated in Fig. 2.c, which describe the bremsstrahlung of additional small- 
X gluons in one step of the evolution. By themselves, such processes are important in the dilute 
regime at relatively low energy (or, for a given energy, at relatively high transverse momenta 
k± ^ Qs), where they generate the n-point correlation functions with n > 2 from the dominant 
2-point function (the gluon distribution). But once generated, the higher n-point functions are 
rapidly amplified by their subsequent BFKL evolution (the faster the larger is n) and eventually 
play an important role in the non-linear dynamics leading to saturation. Thus, such splitting 
processes are in fact important for the evolution towards high gluon density, as first observed in 
numerical simulations^^ of Mueller's 'dipole picture'^'' and recently explained in Refs. 

3 QCD scattering amplitudes from statistical physics 

Equations including both merging and splitting in the limit where the number of colors Nc is 
large have recently became available^, but these are still quite complicated and their exploration 
is only at the beginning (see my next presentation^^). Still, as we shall argue now, the asymptotic 



behaviour of the corresponding solutions — where by 'asymptotic' we mean both the high-energy 
Hmit T ^ oo and the weak coupling hmit ^ — can be a priori deduced from universahty 
considerations relating high-energy QCD to problems in statistical physics^. 

To that aim, it is convenient to rely on the event-by-event description ^ of the scattering 
between the external dipole and the hadronic target, and to use the large-iV^ approximation to 
replace the gluons in the target wavefunction by color dipolcs (which is indeed correct in the 
dilute regime). Then, the scattering amplitude in a given event can be estimated as 



where is the scattering amplitude between two dipoles with comparable sizes and nearby 
impact parameters, and /r(r, b) is the occupation number for target dipoles with size r at impact 
parameter b. Since, clearly, / is a discrete quantity: / = 0, 1, 2, . . . , so is also the scattering 
amplitude in a given event: T is a multiple integer of a^. The estimate (1) is based on the single 
scattering approximation, and thus is valid in the dilute target regime, where T ^ 1. 

In this dipole language, the 2^4 gluon splitting depicted in Fig. 2.c is tantamount to 1 ^ 2 
dipole splitting, and generates fluctuations in the dipole occupation number and hence in the 
scattering amplitude. Thus, the evolution of the amplitude TT{r, b) with increasing r represents 
a stochastic process characterized by an expectation value {T{r, 6)),- ~ {f{r, b))r-i and also by 
fluctuations 5T ~ oil5f ~ a^T, where we have used the fact that 5f ~ ^/f for fluctuations 
in the particle number. These fluctuations are relatively important (in the sense that 5T ^ T) 
only in the very dilute regime where (/) ^ 1, or (T) ^ a^. 

3.1 The mean field approximation (BK, FKPP & geometric scaling) 

Unitarity corrections in the form of multiple scattering start to be important when T ~ 1; 
according to Eq. (1), this happens for dipole occupation numbers of order l/a^. Consider 
first the formal limit 0, in which the maximal occupation number iV ~ l/c^s becomes 

arbitrarily large. Then one can neglect the particle number fluctuations and follow the evolution 
of the scattering amplitude in the mean field approximation (MFA). This is described by the BK 
equation^' a non-linear version of the BFKL equation which, mediating some approximations, 
can be shown ^ to be equivalent to the FKPP equation The latter represents the MFA for the 
reaction-diffusion process A+A and related phenomena in biology, chemistry, astrophysics, 
etc. (see^^ for recent reviews and more references), and reads 



in notations appropriate for the QCD problem at hand: T(/>, r) = {T{r))r and p = ln(rQ/r^), 
with ro a scale introduced by the initial conditions at low energy. Note that weak scattering 
(T ^ 1) corresponds to small dipole sizes (r ^ l/Qs), and thus to large values of p. In 
momentum space, p In/c^. The three terms on the r.h.s. of Eq. (2) describe, respectively, 
diffusion, growth and recombination. Together, the first two terms represent an approximate 
version of the BFKL dynamics, while the latter is the non-linear term which describes multiple 
scattering and ensures that the evolution is consistent with the unitarity bound T < 1. In fact, 
r = 1 is clearly the high-energy limit of the solution to Eq. (2). 

The solution Trip) to Eq. (2) is a front which interpolates between two fixed points : the 
stable fixed point T = 1 (the unitarity limit) at p ^ — oo, and the unstable fixed point T = at 
p ^ oo (see Fig. 3). The position of the front, which marks the transition between strong and 
weak scattering, defines the saturation scale: Psij) = ln(rQQ^(r)). With increasing r, the front 
moves towards larger values of p, as illustrated in Fig. 3. 

The dominant mechanism for front propagation is the BFKL growth in the tail of the 
distribution at large p : the front is pulled by the rapid growth of a small perturbation around 



r,(r,6) ~ alU{r,b), 
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drT{p,T) = dlT{p,T) + T(p,r) - T\p,t) 
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Figure 3; Evolution of the continuum front of the Balitsky-Kovchegov equation with increasing r. 



the unstable state^^. In view of that, the velocity of the front A = dps/dr is fully determined by 
the linearized version of Eq. (2) (i.e., the BFKL equation), which describes the dynamics in the 
tail. By solving the BFKL equation one finds^'''^'^ that, for p > Ps{t) and sufficiently large r. 

Trip) - e'^'^^'e-^" = e-^^P~P^^^)) , p,{t) = ca,r, (3) 

where otg = agNc/TT, 7 = 0.63.., and c = uj/j = 4.88... From Eq. (3) one can immediately 
identify the velocity of the front in the MFA as Aq = ca^. Since Q'^{t) ~ Qq e^""^, it is furthermore 
clear that Aq plays also the role of the saturation exponent (here, in the MFA). 

According to Eq. (3), the scattering amplitude depends only upon the difference p — Ps{t) : 
Tt-{p) = A[p — Ps{t)Y This is an exact property of the FKPP equation (at sufficiently large r) 
and expresses the fact that the corresponding front is a traveling wave which propagates without 
distortion In QCD, this property is valid only within a limited range, namely for^'' 

< p- Ps{t) < Ps{t) (4) 

(the "geometric scaling window"; see below), because of the more complicated non — locality of 
the BFKL, or BK, equations. When translated to the original variables r and r, this property 
implies that the dipole amplitude scales as a function of r'^Q'^^r) : {T{r))T- ^(r^(5g(r)). This 
is the property originally referred to as geometric scaling ^^'^'^'^^ , and which might explain a 
remarkable regularity observed in the small-x data for DIS at HERA. Namely, for x < 0.01, 
the total cross-section (T^*p(x, Q'^) for the absorbtion of the virtual photon shows approximate 
scaling as a function of Q"^ /Ql{x), with Ql{x) oc {1/x)^ and A ~ 0.3 from a fit to the data. This 
measured value of A is quite far away from the above prediction Aq = cus ~ 1 of the BFKL 
equation; but after including the NLO corrections to the BFKL equation (see Ref. for details), 
the ensuing, improved, theoretical prediction decreases indeed to a value close to 0.3. 

3.2 The effects of fluctuations 

What is the validity of the mean field approximation ? We have earlier argued that the gluon 
splitting processes (cf. Fig. 2.c) responsible for dipole number fluctuations should play an 
important role in the dilute regime. This is further supported by the above considerations on 
the pulled nature of the front: Since the propagation of the front is driven by the dynamics in 
its tail where the fluctuations are a priori important, the front properties should be strongly 
sensitive to fluctuations. This is indeed known to be the case for the corresponding problem in 
statistical physics as it can be understood from the following, qualitative, argument: 

Consider a particular realization of the stochastic evolution of the target, and the corre- 
sponding scattering amplitude, which is discrete (in steps of q^). Because of discreteness, the 
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Figure 4: Evolution of the discrete front of a microscopic event with increasing rapidity r. The small blobs are 
meant to represent the elementary quanta of T in a microscopic event. 



microscopic front looks like a histogram and thus is necessarily compact : for any r, there is 
only a finite number of bins in p ahead of ps where T^- is non-zero (see Fig. 4). This property 
has important consequences for the propagation of the front. In the empty bins on the right 
of the tip of the front, the local, BFKL, growth is not possible anymore (as this would require 
a seed). Thus, the only way for the front to progress there is via diffusion., i.e., via radiation 
from the occupied bins at p < ptip (compare in that respect Figs. 3 and 4). But since diffusion 
is less effective than the local growth, we expect the velocity of the microscopic front (i.e., the 
saturation exponent) to be reduced as compared to the respective prediction of the MFA. 

To obtain an estimate for this effect we shall rely on the universality of the dominant 
asymptotic (r ^ oo and N = 3> 1) behaviour, which has been observed in the context 

of statistical physics and justified by Brunet and Derrida through the following, intuitive, 
argument: For a given microscopic front and ^ 1, the MFA should work reasonably well 
everywhere except in the vicinity of the tip of the front, where the occupation number / becomes 
of order one and the linear growth term becomes ineffective. (Note that, in QCD, / ~ 1 



corresponds to T 



~ a: 



which is precisely where one expects the fluctuation effects to become 



important.) Accordingly, Brunet and Derrida suggested a modified version of the FKPP equation 
(2) in which the 'BFKL-like' growth term is switched off when T < a1 : 



drT{p,T) 



al)T{l-T). 



(5) 



By solving this equation in the linear regime, they have obtained the first correction to the front 
velocity as compared to the MFA (in notations adapted to QCD; see Ref. ^ for details): 
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(6) 



where the numbers c ~ 4.88 and k ~ 150 are fully determined by the linear (BFKL) equation. 
In QCD, the same result has been first obtained through a different but related argument by 
Mueller and Shoshi ^. Note the extremely slow convergence of this result towards its mean 
field limit: the corrective term vanishes only logarithmically with decreasing 1/N = a^, rather 
than the power-like suppression usually found for the effects of fluctuations. This is related 
to the high sensitivity of the pulled fronts to fluctuations, as alluded to above. The merely 
logarithmic dependence of Eq. (6) upon the value of the cut-off also explains its universality: a 
renormalization — > Aa^ of the latter does not change the dominant correction in Eq. (6). 

But although it becomes an exact result in QCD in the formal limit — > 0, the estimate in 
Eq. (6) is clearly useless for any practical application, because of the very slow convergence of the 
expansion there. Still, this has the merit to demonstrate that the effects of the fluctuations are 



potentially large, which invites us to critically reexamine the results previously obtained from 
the BK equation. In fact, from the correspondence with statistical physics, we also know^ that 
the geometric scaling property of the BK solution will be eventually washed out by fluctuations 
at sufficiently high energy. But in order to understand how fast this actually happens (i.e., up to 
what energy one should expect geometric scaling to be a good property) , and also to estimate the 
saturation exponent for realistic values of a^, one needs to solve the actual evolution equations 
in QCD at large- A'c ^' to be described in my next contribution to these Proceedings^^. 
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